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ABSTRACT 

The MAGIC collaboration has observed very high energy gamma ray emission from the AGN 1ES 2344+5 14. 
A gamma-ray signal corresponding to an llcr excess and an integral flux of (2.38 ± 0.30 s t a t ± 0.70 sys t) x 
10 -11 cm~ 2 s^ 1 above 200 GeV has been obtained from 23.1 hours of data taking between 2005 August 3 
and 2006 January 1. The data confirm the previously detected gamma-ray emission from this object during 
a flare seen by the Whipple collaboration in 1995 and the evidence (below 5er significance level) from long- 
term observations conducted by the Whipple and HEGRA groups. The MAGIC observations show a relatively 
steep differential photon spectrum that can be described by a power law with a photon index of a = —2.95 ± 
0.12 sta t ± 0.2 syst between 140 GeV and 5.4 TeV The observations reveal a low flux state, about six times 
below the 1995 flare seen by Whipple and comparable with the previous Whipple and HEGRA long-term 
measurements. During the MAGIC observations no significant time variability was observed. 
Subject headings: gamma rays: observations, BL Lacertae objects: individual (1ES 2344+514) 
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All but one of the detected extragalactic very high energy 
(VHE) gamma (7) ray sources so far are active galactic nu- 
clei (AGN) of the BL Lac type. These objects are charac- 
terized by a highly variable electromagnetic emission rang- 
ing from radio to 7-rays, and by continuum spectra domi- 
nated by non-thermal emission that consist of two distinct 
broad components. While the low energy bump is thought 
to arise dominantly from synchrotron emission of electrons, 
the origin of the high-energy bump is still debated. Leptonic 
models ascribe it to inverse Compton processes that either 
up-scatter synchrotron photons (synchrotron-self Compton 
[SSC] models, Marscher & Gear 1985, Maraschi et al. 1992), 
or to external photons that or iginate from the accretion disk 
dDermer & Schlickeiser| [l993), from nearby mass ive stars, or 
are re flected into the jet by surrounding material ( Siko ra et ail 
1994). In hadronic models, interactions of a highly rela- 
tivistic jet outfl ow with ambient matter (iDar & Laor l 119971: 
iBedna rek 1993), proton-induced cascades (Mannheim] 19931) . 
synchrotron radiation off protons (proton synchrotron blazar; 
Aharonian 2000; Miicke & Protheroe 2001), or curvature ra- 
diation, are responsible for the high energy photons. The 
prime scientific interest in BL Lac objects is twofold: (1) to 
understand the VHE 7-ray production mechanisms, assumed 
to be linked to the massive black hole in the center of the 
AGN, and (2) to use the VHE 7-rays as a probe of the extra- 
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galactic background light (EBL) spectrum between about 0.3 
to 30 /Ltm wavelength. In the past, most of the VHE 7-ray 
emitting AGN were discovered during phases of high activity, 
biasing our current observational record towards high emis- 
sion states. Although these sources also show variability in 
the X-ray, optical, and radio domain, the VHE variability is 
observed to often be the most intense and violent one. While 
fast variability on the timescale of 10 min utes has been ob - 
served for Mkn 501 in the X-ray domain dXue & Cuill2005l) . 
flux doubling times well bel ow 5 minutes were recently also 
found in the VHE domain dGaidos etalj|1996t lAlbert et ail 
l2007cl) . Many of the observed AGN are presumably visi- 
ble only during a state of high activity. It still remains an 
open question whether these sources are only temporarily ac- 
tive and are completely inactive between times of flaring, or if 
there also exists a state of low but continuous 7-ray emission. 
In addition, the temporal and spectral properties of such a low 
VHE 7-ray emission state is mostly elusive as of to date. It is 
quite conceivable that, compared to a low state, the flare emis- 
sion state is either due to a different population of accelerated 
particles or originates from a different region in the AGN, or 
both. 

In the first year of operation of the Major Atmospheric 
Gamma Imaging Cerenkov (MAGIC) Telescope a program 
has been started to search for new low and medium redshift 
blazars emitting at VHE 7-rays (Albert et al. 2006b; 2006c; 
2007a). In addition, known VHE AGNs were monitored in 
order to study common features of their 7-ray emission, as 
well as the properties of the low-emission state (Albert et al. 
2006a; 2007b; 2007c). 

A good candidate for detailed studies is 1ES 2344+514. 
This AGN belongs to a type of blazars in which the syn- 
chrotron emission peaks at UV/X-ray frequencies (the so- 
called high-energy peak BL Lacs [HBLs] e.g. Urry & 
Padovani 1995), as opposed to the blazars with the syn- 
chrotron peak located at IR/visible frequencies. Along with 
Mkn 501 and H 1426+428, it represents extreme BL Lac ob- 
jects, in which the synchrotron pea k energy exceeds 10 keV , 
in particular during strong flares (Costam ante et al.l 2001). 
1ES 2344+5 14 was d etected during the Einstein Slew Sur- 
vey (lElvis et alJll992l) in the energy range be tween 0.2 and 
4 keV It was identified as a BL Lac object by Perlm arTet alj 
( 1996), who also determined a redshift of z = 0.0 44. Its black 
hole m ass was estimated to be 10^ 8 - 80±0 16 )M g „n (iBarfh et all 
12003b . Early BeppoSAX observations dGiommi etalj|2000t) 
revealed a large 0.1 — 10 keV flux variability on timescales 
of a few hours. Follow-up observations in 1998 found the 
object in a very low state with the synchrotron peak shifted 
by a factor of 20 towards lower energies and t he correspond- 
ing int egral flux decreased by a factor of 4.5. Giom mi et alj 
(2000) interpreted the observations with one electron popula- 
tion being responsible for the steady low energy synchrotron 
emission and another electron component producing higher 
energy X-rays with high time variability. The latter compo- 
nent should be responsible for VHE 7-ray emission via in- 
verse Compton (IC) scattering. EGRET did not detect any 
signal from 1ES 2344+514, giving an upper limit of 3.4 x 
10~ n erg cm" 2 s -1 at its peak response energy of 300 MeV 
dFichtel et al.lll994l) . During the winter of 1995/1996, the 
Whipple collaboration reported a 5.8er excess signal from 
1ES 2 344+514 above 350 G eV from 20.5 hours observation 
time dCatanese et al.l [19981) . The observed flux was highly 
variable, with the most significant signal occurring during 
a flare on 1995 December 20, while all the remaining data 



combined led to an only marginal (4<t) excess, i.e. below the 
canonical detection limit used in ground-based VHE 7-ray as- 
tronomy. The 0.8 — 12.6 TeV differential spectrum measured 
by the Whipple collaboration during the fl are had a power- 
law index of -2.54 ± 0.17 sta t ± 0.07 syst dSchroedter et alJ 
2005). One year later another search did not reveal any VHE 
7-ray emission. The HEGRA collaboration also searched 
for VHE 7-ray emission above 800 GeV. A deep exposure 
of 72.5 h indicated a signal at a significance level of 4.4<r 
dAharonian et alj|2004l) . 

Here we present MAGIC telescope observations of 
1ES 2344+514. We briefly discuss the observational tech- 
nique used and the implemented data analysis procedure, de- 
rive a VHE 7-ray spectrum of the source, and put the results 
into perspective with other VHE 7-ray observations of this 
AGN. An SSC model is used to describe the wide-range spec- 
tral energy distribution (SED). 

2. OBSERVATIONS 

The observations were performed between 2005 August 3 
and 2005 September 29, and between 2005 November 11 
and 2006 January 1, using the MAGIC Telescope on the Ca- 
nary island of La Palma (28.8°N, 17.8°W, 2200 m above sea 
level), from where 1ES 2344+514 can be observed at zenith 
distances above 24°. The essential parameters of the currently 
largest air Cherenkov telescope are a 17 m segmented mirror 
of parabolic shape, an f / D of 1.05 and a hexagonally shaped 
camera with a field of view (FOV) of ks 3.5° mean diameter. 
The camera comprises 576 pixels composed of hemispherical, 
six dynode photomultipliers augmented in sens itivity by a dif- 
fuse l acquer doped with a wavelength shifter dPaneque et alj 
2004) and by so-called light catchers. In separate measure- 
ments a total gain of 2 has been determined. 180 pixels of 
0.2° surround the inner section of the camera, which con- 
sists of 394 pixels of 0.1° (= 2.2° FOV). The trigger 
is formed by a coincidence of > 4 neighboring pixels. The 
overall Cherenkov photon (300 — 650 nm) to photoelectron 
conversion ratio is 0.15 ± 0.02. The point spread function 
(PSF) of the main mirror is a « 0.04°, while 90% of the light 
of a source at infinity is focussed onto a disk with 0.1° 0. Fur- 
ther details of the telescope parameters and performan ce can 
be found in lBaixeras et all d2004l) : ICortina et all d2005l) . 

1ES 2344+514 was observed for 32 hours in total, dis- 
tributed over 27 days between 2005 August and the first days 
of 2006 January at zenith angles ranging from 23 ° to 38° . The 
obser vations were carried out in wobble mode dFomin et all 
1994), i.e. by alternatingly tracking two positions at 0.4° off- 
set from the camera center. This observation mode allows a 
reliable background estimation for point sources. 

Simultaneous i?-band observations of 1ES 2344+514 were 
conducted in the framework of the Tuorla Observatory blazar 
monitoring program 23 with the KVA 35 cm telescope 24 on La 
Palma and the 1.03 m telescope at Tuorla Observatory, Fin- 
land. 

3. DATA ANALYSIS 

The data analysis was carried out using the stan- 
dard MAGIC analys is and reconstruction software 
dBretz & Wagner! 120031) . After calibration dGaug et all 
120051) . the images were cleaned by requiring a minimum 
number of seven photoelectrons (core pixels) and fi ve pho- 
toelectrons (boundary pixels), see e.g. iFeganl (1 19971) . These 

23 See http://users.utu.fi/kani/lm/ 

24 See http://tur3.tur.iac.es/ 
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tail cuts are scaled accordingly for the larger size of the 
outer pixels of the MAGIC camera. The data were filtered 
by rejecting trivial background events, such as accidental 
noise triggers, triggers from nearby muons or data taken 
during adverse conditions (low atmospheric transmission, 
car light flashes etc.). Light clusters, either from large angle 
shower particles or from the night sky light background 
(stars), well separated from the main image, were removed 
from the images. For the events included in the analysis, 
the mean trigger rate was required to be constant within 
?» 20%. In order to improve the comparability of the two 
data sets from summer and winter 2005 we restricted the 
maximum zenith angle to < 34°. From the remaining events, 
corresponding to 23.1 h obser vation time, image parameters 
were calculated dHillasI [T985) such as WIDTH, LENGTH, 
SIZE, CONC, and M3LONG, the third moment of the light 
distribution along the major image axis. For the 7/hadron 
separation a multidimensional classification procedure based 
on the random fo rest method was employed (Breimar feOOU 
iBock et ail 120041) . The separation procedure was trained 
using a s ample of Monte Carlo (MC) generated 7-ra y shower 
images dHeck et all 119981; iMajumdar etail 120051) on the 
one hand and about 1% randomly selected events from the 
measured wobble data representing the hadronic background 
on the other hand. The MC 7-ray showers were generated 
between zenith angles of 24° and 34° with energies between 
10 GeV and 30 TeV. Every event was assigned a parameter 
called hadronness (h), which is a measure for the probability 
that it is a hadronic (background) event. The final separation 
was achieved by a cut in h. The same cut procedure was 
applied to the final 1ES 2344+514 sample. The arrival direc- 
tions of the showers in equ atorial coordinates were calculated 
using the DISP method (jFomin et al.1 Il994t iLessard etail 
1200 it iDomingo-Santamarfa et al.ll200l . The energy of the 
primary 7-ray was reconstructed from the image parameters 
again using the random forest method and taking into account 
the full instrumental energy resolution. 
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FIG. 1.— Sky map for 1ES2344+514 produced with a DISP analysis: 
The figure shows the (background-subtracted, see e.g. Rowell 2003) excess 
events above 300 photoelectrons, corresponding to a 7-ray energy of as 180 
GeV. The sky map has been smoothed using a two dimensional Gaussian with 
cr = 0.1°, roughly corresponding to the 7 PSF of the MAGIC telescope for 
point sources (indicated by the white circle). The colors encode the number 
of excess events in units of 10 — 5 sr _1 . The black cross marks the expected 
source position. 
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FIG. 2. — 9 2 plot for the 1ES 2344+514 observations. On-source events are 
given as black symbols, while red symbols represent off-source background. 
A cut at 9 2 < 0.04 deg 2 selects 528 7 events at a significance level of llcr. 
The plot has been prepared for SIZE > 300 photoelectrons, corresponding 
to « 180 GeV. 

Fig. Q] shows a w 1.8° x 1.8° section of the sky around 
the 1ES 2344+514 position. The nominal source position is 
marked by a cross. A clear excess is visible in the data, 
the maximum of which is located at (RA,Dec)=(23 h 46 m ± 
0™4, 51°42.'6 ± 1.'2) (the errors only include the determina- 
tion accuracy of the position). The extension of the excess and 
the small deviation from the nominal position are consistent 
with the PSF and the tracking error of w 1™5 of the telescope 
dBretz et al. 2003), respectively. To calculate the significance 
of the observed 7-ray excess, the squared angular distance 9 2 
between the reconstructed shower direction and the object po- 
sition (0 2 = 0) as shown in Fig.|2]is used. In this represen- 
tation, the background is expected to be flat for the case of a 
very large diameter camera. In the analysis, three background 
regions of the same size chosen symmetrically with the source 
position around the camera center were used for a simultane- 
ous determination of the background. The background control 
data sample was normalized to the on-source sample between 
0.12 deg 2 < 8 2 < 0.24 deg 2 . The reason of the slow but 
steady drop in the background is the drop in acceptance to- 
wards the camera boundary. The observed excess signal of 
528 events below 2 < 0. 04 deg 2 correspo nds to an llcr ex- 
cess according to eq. 17 in Li & Ma] (Il983h . An independent 
analysis using other cuts, a different reconstruction algorithm 
and a different 7/hadron optimization procedure, revealed a 
comparable (within statistics) significance. While for the sky 
map and the 2 plot a fixed, tight h cut was applied, the fi- 
nal separation for the spectral analysis and the light curve was 
done using a looser, energy-dependent cut in h, requiring that 
about 60% of the MC 7 events survive. 

As the analyzed data comprise 21 observation nights, it 
is possible to check the light curve for possible flux vari- 
ability. On a diurnal basis, the > 200 GeV light curve 
(Fig. [3} shows small changes and trends beyond those 
expected from statistical fluctuations. The structure ob- 
served during MJD 53580—53600 is compatible with a 
constant-flux ansatz (x 2 /dof= 6.1/6), while from MJD 
53726.82-53726.90 a flux of 2.4ct above the average flux in- 
ferred from the surrounding days MJD 53720-53740, (1.8 ± 
0.6) • 10~ 11 cm^ 2 s" 1 (x 2 /dof=4.9/7), was found. Note that 
the probability for finding such an excess in the 21 observa- 
tion nights is around 34.4%. No significant variability within 
this observation night, encompassing 1.13 hours of effective 
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FIG. 3. — VHE (E > 200GeV) light curve (upper panel) and simulta- 
neous optical (/2-band) light curve for 1ES 2344+514. The dashed line in 
the VHE light curve indicates the average flux level of (2.38 ± 0.30) X 
10~ 11 cm _2 s- 1 (x 2 od = 21.2/20). Note that the contribution of the 
host galaxy to the o ptical brightness is non-negligible and given as 3.7 mjy 
Wilsson et aIT2 006). 
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Obs. time 
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X7"df 


(MJD) 




(lO" 11 cm- 2 s- 1 ) 




53585 - 53597 


6.37 h 


3.02 ± 0.50 


6.1/6 


53610 - 53642 


8.06 h 


1.87 ±0.52 


0.4/3 


53700 - 53736 


8.66 h 


2.20 ±0.51 


12.0/9 


Combined 


23.09 h 


2.38 ± 0.30 


21.2/20 



TABLE 1 

Integral fluxes above 200 GeV in the individual observation 
periods and reduced x 2 /dof of a fit with a constant-flux 
ansatz in the respective observation periods. the given 
errors are statistical errors only. 

observation time, can be claimed. 

The observation time can be split into three observation pe- 
riods (Table[T|i. Together with the VHE 7-ray light curve, an 
i?-band optical light curve is shown. Simultaneou s X-ray data 
are o nly available from the ASM instrument 25 dLevine et alj 
1996) on board the Rossi X-ray Timing Explorer, the sen- 
sitivity of which, however, would be hardly sufficient to 
resolve the expected 2 — 10 keV flux even during flaring 
states of 1ES 2344+5 14, like those observed with BeppoSAX 
dGiommi et al.ll2000h . 

Summing up all the data we determined an integral flux 
above 200 GeV of 

F(E > 200GeV) = (2.38±0.30 st at±0.70 syst )xl0" 11 cm- 2 s 

The relatively large systematic error is a consequence of the 
steep spectral slope (see below). The main contributions to 
the systematic error are the uncertainties in the atmospheric 
transmission, the reflectivity (including stray-light losses) of 
the main mirror and the light catchers, the photon to photo- 
electron conversion calibration, and the photoelectron collec- 
tion efficiency in the photomultiplier front end. Also, MC 
uncertainties and systematic errors from the analysis methods 
contribute significantly to the error. The above quoted flux 
corresponds to (10 ± 1)% of the integral Crab Nebula flux in 
the same energy range. 

During our observations we also checked the optical vari- 
ability. When cor recting for the con tribution of the host 
galaxy of 3.7 mjy (Nilsson et al. 2006), variations in the op- 
tical light curve around the average brightness of w 15% are 
seen, which are significant given the small errors (< 5%) of 
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TABLE 2 

Differential flux of 1ES 2344+5 14 along with statistical and 

SYSTEMATICAL ERRORS. 

the data points. Possible VHE 7-ray variations on a compa- 
rable level are below the sensitivity of MAGIC on the given 
timescale. 

For each of the three observation periods photon spectra 
were determined. These are well described by simple power 
laws between 140 GeV and at least 1.0 TeV and are, within 
errors, compatible with no change in the spectral index. Fi- 
nally, all data were combined for the calculation of a differ- 
ential photon spectrum (Table [2]). The reconstructed spec- 
trum after unfolding w i th the instrumental ene rgy resolution 
(lAnykeyev et al J fl99 1 : Mizobuchi et alJ 120051) is shown in 
Figure [4] A simple power law fit to the data between 140 
GeV and 5.4 TeV yields 
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with a reduced x 2 /dof of 8.56/5, indicating a reasonable de- 
scription of the data by the fit. For comparison, the Whip- 
ple measurement of the 1ES 2344+5 14 spectrum du ring the 
flare of 1995 Decembe r 20 (ISchroedteret al.ll2005l) and the 
Crab Nebula spectrum (Wagner et al. 2005) as obtained with 
MAGIC are also shown in Fig. [4] Note that the integral flux 
F(E > 970 GeV) = (0.82 ± 0.09) x 10~ 12 cm" 2 s" 1 is in 
very good agre ement with the HEG RA measurements from 
1998 to 2002 (lAharonian et al.ll2004l) . 
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FIG. 4. — Differential photon spectrum for 1ES 2344+514 as measured 
with MAGIC (solid curve). The gray band represents systematic errors com- 
ing from varying the 7 efficiency in the determination of the spectrum. The 
Crab nebula spectrum as measured with MAGIC is also shown (gray dotted 
curve; small gray numbers indicate the fraction of Crab nebula flux for the 
1ES 2344+ 514 flux points) and the Whipple flare spectrum (dashed curve) as 
reported by Schroedter etal] <2005h . 



4. DISCUSSION 
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FIG. 5.— Overall VHE light curve for 1ES 2344+514. Data for 1996-2002 were collected from ICatanese et alj [\99H> : iBadran & Weekes! HMHr . 
lAharoni an et al. ( 2004). The 1995 December 20 flare has been excluded for clarity. Due to the different energy thresholds of the included observations, integral 
fluxes > 350 GeV were considered and the data were normalized to Crab flux units, extrapolating data points with higher flux levels assuming 1ES 2344+514 
was observed in a quiescent state; the spectral behavior found in this paper was used for this purpose. 

We observed a highly significant VHE 7-ray emission from 
the blazar 1ES 2344+514. The flux exhibited no significant 
variations on a timescale of days, with one night showing a 
higher flux level (by a factor of 2) as compared to the sur- 
rounding nights; such a variation is statistically expected to 
occur with a 34.4% probability in 21 observation nights. The 
observed flux was lower by a factor of k 6 than the one ob- 
served by the Whipple collaboration during a flare on 1995 
December 20. The spectrum is softer than the Crab Nebula 
spectrum and also softer than the flare spectrum observed by 
the Whipple collaboration. 



4.1. Long-term VHE Light Curve 

In Fig. [5] we show a light curve including all reported VHE 
7-ray measurements and upper limits for 1ES 2344+514. 
These data have been normalized to an integral Crab flux 
F{E > 350 GeV); the fluxes given in the literature were 
extrapolated, if necessary, using the spectral index found in 
this paper (in the following, "Crab units" refer to this energy 
threshold). Except for the 1995 December 20 flare and the 
MAGIC data, all these measurements are on the sensitivity 
level of the respective instruments. Therefore, none of the lat- 
ter data points exceeds a significance of 4.3rj. All reported 
observations with significances below 2. Oct («95% probabil- 
ity) were converted to 99% upper flux limits. 

In 1995/1996, Whipple discovered 1ES 2344+514 at a flux 
level of (0.11 ± 0.05) Crab units at E > 350 GeV, except 
for the 1995 D ecember flare, when (0.63 ±0.15) Crab units 
were obtained (ICatanese et al.lfT998h . Follow-up observations 
by Whipple and HEGRA in 1996-1998 yielded upper lim- 
its of 0.08 and 0.12 Crab units, respectively. In 1998 and 
2002, the object was observed for almost 60 h by HEGRA re- 
sulting, when combined, i n a flux of ( 0.042 ± 0.012) Crab 
units at E > 930 GeV dTluczvkontl 120031) . which trans- 
lates to (0.053 ± 0.015) Crab units when extrapolating to 
E > 350 GeV From observations of 1ES 2344+514 in 2002, 
the Whipple group could infer a low flux le vel of < 03 Crab 
units with a marginal significance of 3.1ct dBadran & Weekes] 
12001 at E > 400 GeV 

While the Whipple and HEGRA measurements allowed 
to conclude on a emission level of < 11% Crab units only 
after long observation times, the MAGIC observations ob- 
tained in this paper are the first time-resolved measurements 
at this emission level for 1ES 2344+514. We find the flux 
of 1ES 2344+514 to be (0.054 ± 0.006) Crab units for E > 
350 GeV, which is well in line with the HEGRA 1997-2002 
evidence. 

In previous observations of 1ES 2344+514 it was not pos- 



sible to infer temporal characteristics of the found VHE 7- 
ray emission level. With MAGIC, this level can be detected 
with only a few hours of observations, enabling studies of the 
VHE 7-ray variability properties of this object over a sig- 
nificant part of its dynamical range. Thus, 1ES 2344+5 14 
adds to the small group of blazars for which such studies 
are no w possible o n a diurnal basis — M kn 421 (lAlbert et alJ 
120063). Mkn 501 J Albert et alJ 120073). and PKS 2155-304 
( lAharonian et alj|2005h ~ 

4.2. Intrinsic Energy Spectrum 

Having to traverse a cosmological distance correspond- 
ing to a redshift of z = 0.044, the 7-rays emitted by 
1ES 2344+514 interact with the low energy photons of the 
EBL (see, e.g., Nikishov 1962; Gould & Schreder 1966; 
Hauser & Dwek 2001). The predominant reaction 7vhe + 
7ebl — ► e + e~ leads to an attenuation of the intrinsic spec- 
trum dN/ dEiatr that can be described by 

dN/dE ohs = dN/dE intI • exp[-r 77 (i;, z)] (1) 

with the observed spectrum diV/di? bs, and the energy- 
dependent optical dept h t^(E, z). Here we use the 
"best-fit" model of Kne iske et alJ (120041) . which yields an 
EBL spectrum tha t agrees with alternative models, e.g. 
IStecker et alJ ((2006). Using this EBL spectrum and a state- 
of-the-art cosmology (flat universe, Hubble constant Hq = 
72 kins -1 Mpc -1 , matter density fl m = 0.3, and dark energy 
density Q\ = 0.7), we calculate the optical depth t 77 for the 
distance of 1ES 2344+514. Thereby we use the numeri cal in- 
tegration given by eq. 2 in iDwek & Krennri ch (2005). The 
reconstructed intrinsic source spectrum is shown along with 
the measured spectrum in Fig. [6] The intrinsic source spec- 
trum can be described by a simple power law of the form 



dN 

dE intr 



(2.1±1.2 stat ±0.5 syst )-10- n E 



-2.66±0. 



TcV cm 2 s 



500 GeV 



between 140 GeV and 5.4 TeV (x\ oi = 0.68/5). The spec- 
trum shows a tendency to flatten towards low energies. A 
fit with a logarithmic curvature term, which corresponds to a 
parabolic law in a \og(E 2 dN / dE) vs. \ogE representation 
(power law index a — > a + 2b\og(E / E a ); Massaro et al. 
2004), shows a clear curvature and enables locating the peak 
at £ P cak = E a ■ I0( 2 " a )/( 2b ) = (202 ± 174) GeV. The peak is 
obviously badly determined as the turnover of the spectrum, 
presumably around 200 GeV, is not observed unambiguously. 

4.3. Spectral Energy Distribution 
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FIG. 6. — Measured and intrinsic differential photon spectrum of 
1ES 2344+5 14. Simple power law fits are given by the dot-dashed and dashed 
lines, while the fit with a logarithmic curvature term is given by the dotted 
curve. 

The spectral energy distribution for 1ES 2344+5 14 is shown 
in Fig. [7] Apart from the VHE 7-ray observations by MAGIC, 
Whipple, and HEGRA, the X-ray measurements performed 
by Einstein, ROSAT, and BeppoSAX and an upper limit from 
EGRET are the most relevant measurements for modeling 
the SED. The figure also includes optical and radio measure- 
ments. The latter however cannot be describ ed by the homo- 
geneo us one-zone SSC model provided bv lKrawczvnski et alj 
(2004) that was used here, as the radio emission is thought to 
arise from a larger volume in the jet than the VHE 7-ray emis- 
sion and as self-absorption effects are not accounted for. 

Two BeppoSAX data sets are shown: one was obtained si- 
multaneously with Whipple observations on 1996 December 
5. The second set was taken during a quiescent period of 
1ES 2344+514, and represents a rather low X-ray emission 
state of the object, although it might not necessarily corre- 
spond to its state during the MAGIC measurements. It should 
be emphasized that most of the data points shown in the SED 
were not measured simultaneously, which makes the SED 
modeling very difficult. The input parameters of the homo- 
geneous SSC model are the radius of the spherically assumed 
emission region R, the Doppler factor <5, the magnetic field 
strength B in the acceleration region, the density p of the elec- 
trons responsible for the 7-ray emission as well as two spec- 
tral slopes and a spectral cutoff of the electron spectrum. The 
size of the emission region R was chosen such that it can ac- 
count for day-scale variability along with the Doppler factor 8 



chosen, R < <5_i x t d x 2.48 x 10 16 cm with £_i = 5/10 and 
t c \ in units of days. The parameters used here are specified in 
Table[3] 

The two models shown are to represent both the flare state 
of 1995 and the MAGIC observations in 2005. They differ 
in the following model input parameters: (1) in the Doppler 
factor, (2) in the magnetic field strength, and (3) the minimum 
and maximum electron energy. In addition, the form of the 
electron spectrum differs for the two flux states. The radius of 
the emission region was kept constant at 10 14 m. While most 
of the model parameters are compatible with the parameter 
space spanned by other models for blazars (e.g. Kino et al. 
2002; Giommi et al. 2002), the magnetic field strength found 
here is rather low. 

In conclusion, we note that the presented SED models 
are rather speculative, given the non-simultaneity of the cur- 
rently available data. Future multiwavelength campaigns on 
1ES 2344+5 14, exploiting the enhanced sensitivity of the new 
im aging air Cherenkov telescope installations, will hopef ully 

SSC model input Low state Flare SED 



parameter SED (MAGIC) 



Doppler factor <5 


8.4 


15.2 


Magnetic field strength B 


0.095 G 


0.075 G 


Emission region radius R 


10 16 cm 


10 16 cm 


Electron density p 


0.025 erg cm -3 


0.025 erg cm" 3 


J B min [log£(eV)] 


9.1 


8.9 


£ max [logE(eV)] 


11.6 


11.9 


E br [logE(eV)] 


10.9 


10.9 


n(E min < E < E hr ) 


-2.2 


-2.2 


n(E hT <E< £ max ) 


-3.2 


-3.2 



TABLE 3 

M odel input parameters use d with the SSC model provided by 
IKrawczynski et al] ( 120041) for describing the MAGIC low 

EMISSION STATE AND THE 1995 WHIPPLE FLARE EMISSION STATE AS 
DEPICTED IN FIG. [7J 



improve this situation. 
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